arXiv: 1505.07443v2 [hep-ph] 3 Aug 2015 


Analytical two-loop soft mass terms of sfermions in Extended GMSB models 


☆ 


Tomasz Jelinski*, Janusz Gluza 

Institute of Physics, University of Silesia, Uniwersytecka f, PL-fO-007 Katowice, Poland 


Abstract 

Analytical two-loop contributions to soft masses of sfermions are derived in the Extended GMSB (EGMSB) model with 
one superpotential coupling between matter and messenger superfields. Analytical results allows to study in detail the 
whole range of the ratio F/M^ of the spurion F-term, F, to the messenger scale M. It is shown that if F/M'^ is of the 
order of 1, then one- and two-loop contributions to soft masses are of the same magnitude and their interplay leads to 
novel sfermion mass patterns. 
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1. Introduction 


Supersymmetric models which involve GMSB mecha¬ 
nism yy] are popular frameworks for analyzing phenomenol¬ 
ogy of BSM scenarios mainly due to their flavour univer¬ 
sality of soft SUSY breaking terms. 

In this letter, we consider in more details the issue of 
sfermion masses in extended GMSB scenario. Let us start 
by reviewing situation within pure GMSB scheme. 

It is well-known that in such s^ups, I-loop A-terms 
and I-loop soft masses of sfermions 4>i are zero at the mes¬ 
senger scale M. The first non-trivial contribution to soft 
masses arise at 2-loop level and can be written as 1,1: 

= CJiMxr/M ( 1 ) 


where gr are gauge couplings for 51/(3), SU{2) and U{1) 
groups, respectively. N is twice the sum of the Dynkin in¬ 
dices of the messenger fields, is the quadratic Gasimir 
related to (fti, while x = F/M'^ is a dimensionless param- 
etei0 which sets the relation between F-term of the spu¬ 
rion, F, and the messenger scale M. Einally, the function 
fgAx), 
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only mildly depends on x such that the overall scale of © 
is set by gf{Mx)'^j (IGtt^)^. This quantity is crucial for the 
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phenomenology of GMSB models. To arrange masses of 
sfermions at the level of a few TeV, such that sparticles 
are accessible at the LHG2 or at one of its successors, and 
to avoid fine-tuning of parameters, the scale Mx has to 
be of the order of 10® GeV. This allows us to consider the 
following relation betweerd M and x 

Mx = X 10® GeV, (3) 

where cj is a numerical constant of the order of 1. Its 
precise value depends on the details of the model (e.g. 
tan/3, /i-term). 

However, the GMSB models are not flexible enough to 
fulfil all phenomenological restrictions and simultaneously 
retain Eor example, it is hard to accommodate for 
125 GeV Higgs boson mass when sfermions are relatively 
light 0. It turns out that the above mentioned issue can 
be alleviated within so-called Extended GMSB (EGMSB) 
models in which one allows for superpotential coupling be¬ 
tween messenger and MSSM fields. Beside addressing en¬ 
hancement of the Higgs boson mass, EGMSB models also 
provide reasonable framework for exploring non-ininimal 
flavour violation scenarios within various GUT models 
i i E [H [13, [H- As a result of the presence of addi¬ 
tional contributions to soft terms, EGMSB models exhibit 
rich and interesting phenomenology in which e.g. light 
stops can be natur ally realized ilillil. 

It is known [13 [l3 that in EGMSB models both I- and 
2-loop soft masses are generated by the messenger-MSSM 
superpotential couplings. The I-loop masses are known 
analytically for any 0 < x < I and one can show that 
they are negative and oc h?x‘^{xMY for a; 1. On 
the other hand only the leading term in a; = F/M^ of the 
2-loop contribution to soft masses is known. It has been 


:^The mass scale Mx is very often written as FfM, see e.g. Q 
and sometimes denoted by see mu. 
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extracted with the h elp of the wave-function renormaliza¬ 
tion method d, [H, [3, 0, 0, EH, which captures only 
contributions of the order of One can show that 

for X <C 1 the 2-loop contributions to soft masses are posi¬ 
tive and oc /{lOnr'^Y. Naively, taking into account 

above considerations one expects that 1-loop contributions 
could start to dominate for x > ~ 10“^ — 10“^. 

Though, such rough estimation does not take into account 
higher powers of x in the expansion of 2-loop masses, but 
also does not encounter numerical factort0 which originate 
from contracting tensors related to the representations of 
chiral superfields 17|, [21 1. 

Let us recall that in many GUT models the messenger 
scale is chosen rather largeo e.g. M = 10^'^ GeV, what 
due to m forces X to be of the order of 10 As men¬ 
tioned above for such choice of x the full form of the loop 
function @ is not necessary - only the leading term is rele¬ 
vant. However, there are also interesting setups, so-called 
Low-Scale Gauge Mediation models with relatively light 
messengers M ~ 10® — 10® GeV which provide interesting 
phenomenology, see e.g. (^ . In this situation x can be 
much larger, even of order 1. In this regime wave-function 
method is not adequate. To provide more accurate es¬ 
timation of soft masses, one needs full analytical 2-loop 
expression. 

Our aim is to calculate such full analytical formula 
for 2-loop soft masses of sfermions which is valid for all 
0 < X < 1, examine its relation to 1-loop soft masses and 
check for which values of x full analytical expressions are 
important. To derive them, we adopt a general method 
presented in 0, 261, which is suitable for 2-loop SUSY 


calculations, and conventions therein. 


2. EGMSB model with one messenger-matter cou¬ 
pling 

We consider a EGMSB-type model in which a pair of 
messengers, Y, Y couple to MSSM fields (/>i, (j )2 via su¬ 
perpotential term hY(j}i(f >2 and SUSY is broken by vev of 
F-term of spurion superfield X. In other words, the rele¬ 
vant part of the superpotential is of the following form: 

Wh = XYY + h<l>i(l)2Y. (4) 

The model has symmetry (j)i -o- (j )2 and, moreover, one 
can assign global t/(l) charges to (jii and Y, Y such that 
those fields cannot mix with each other even after SUSY 
breaking^ An example of such assignment is the following: 


^These numerical factors also enter 1-loop anomalous dimensions 
of chiral superfields Sam. 

the GMSB models gravitino mass is m^j2 ~ F/Mpi = 
xM'^fMp^ where Mp ~ 10^® GeV is the Planck scale. If ((3ll is 
satisfied then 771^/2 10“® GeY/x. Hence in the models with x 1 

the gravitino mass is of the order of 10 eV. Such li ght gravitino may 
be of some importance for the cosmology, see e.g. . 

^For the discussion of the role of such charges in the GUT models 
see e.g. [i,[ll[2||. 


> 0, 901 ^ 902, 901 -f 902 = -qy = 9y and qx = 0. 
Due to the above-mentioned symmetries: (i) both fields 
</>i ,2 obtain the same soft mass (ii) there is no mixing 
mass term (/)i^ 2 +h.c. in the scalar potential, and (hi) there 
is no mass term of the form + h.c. Finally, (iv) one 
can get rid of any complex phase oi h = just by a 

redefinition of messenger fields: Y —?► Y —Ye®’^'*. 

Hence, without loss of generality, we can assume that h 
is real. An example of superpotential interaction (P|) is a 
coupling between up Higgs 4>i = iL„, up quark 4>2 = U 
and messenger Y = Yg in the representation (3, 2)i/e of 
SU{3) X 5't/(2) X C/(l). 

When spurion field gets vev {X) = M + 0^F, all com¬ 
ponents of messenger multiplets obtain masses M, while 
scalar components get extra SUSY breaking masses FYY+ 
h.c. After unitary rotation of the scalar components of 
messengers (Y, Y) to the basis (Y+, Y_) in which their mass 
terms are diagonal, the scalar potential is of the form: 


V = Ml\Y+\^ + Ml\Y_\ 


—=hMcj)i(j) 2 {Y+ -I- Y_) -I- h.c. 

v2 


F^h\\cj^,\^ + \cl^2?)\Yy-Y.\\ 


( 5 ) 


where M|. = M^(l ± x) and x = F/M'^. Let us remark 
that X has an upper bound, related to the upper bound on 
the mass of Y_, i.e. on M_ = My/1 — x. Because Y_ can 
have the same quantum charges as quark superfields hence 
LHC bound on squarks masses ttilhc ~ 2 TeV 0, E3 
apply to Y_ as well. However, for the simplicity, in this 
letter we assume a little more conservative bound: 


A/_ > 10^ GeV, 


( 6 ) 


such that Y_ is significantly heavier than all sfermions. 
Otherwise Y_ may, in some cases, alter stability of the 
scalar potential, what in turn can be dangerous for EWSB. 
Using definition of M_, the condition ([6]) can be rewritten 
in the following form: 


X < 1 — 


/10®‘ GeV 
V M 


2 


( 7 ) 


It ensures that below the scale of 10 TeV one gets only 
MSSM fields. For example, when M = 10® GeV, x has 
to be smaller than about 0.99, but for M = 10^^ GeV, 
that limit is practically equal to 1 and may indicate a 
kind of severe hne-tuning of parameters: x < 1 — 10“^®. 
Additionally, taking into account (I3|) the bound ([7]) can be 
rewritten as: 



( 8 ) 


Let us summarize our discussion. For x <C 1, the 2- 
loop soft masses can be obtained with the help of the 
wave-function renormalization method. We also know that 
the regime x < 1 is interesting for the phenomenology of 
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Low-Scale Gauge Mediation models. But for a; ^ 1 the 
wave-function method is not adequate. Hence, to fill that 
gap, we shall discuss full analytical result for 2-loop masses 
generated by O- 


3. 1-loop soft-terms 


For the completeness, we shall present 1-loop H-terms 
and 1-loop soft masses in the discussed model. We omit 
soft masses of gauginos because at the 1-loop level thw 
are of the same form as in the standard GMSB models [l| . 
The 1-loop soft mass of 4>i is 


m 


2 _ 
h,i — 


Jf_ 

IGtt^ 





(9) 


where the function fh,iix) is dehned as follows: 

fh.ii^) = [(a: - 2) ln(l - a;) - (x-I-2) ln(l-I-a:)] 

= 1-1- -I- -^x"' -I-- (for X <C 1) (10) 

5 14 

On the other hand the H-term is: 

^ ( 11 ) 


where 

A,a(x) 


1 , 1 -I- X 

— In- 

2x \ — X 

1 -I- \x^ + ^x^ - (for X < 1) 
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( 12 ) 


It is clear that both m\ ^ and A diverge when x —>■ 1. The 
reason for that is the massless scalar IT which is present 
in the 1-loop diagrams [3l| and leads to the divergence 
of the following contribution (M^/i^/IGtt^) ln(M^/M^) to 
the self-energy H^^^ (0) at zero momentum s = = 0. To 

ensure that the perturbative expansion in /i^/IGtt^ does 
not break down, the following conditions has to be satisfied 


|1 


x| exp 



(13) 


the pole mass of 4>i. One can show that 1- and 2-loop pole 
masses are of the following formal 


'-h,! 


1 

IGtt^ 


nW(o), 


(14) 


''h,2 


1 

(167r2)2 


n(i)'(o)n(i)(o) -f n(2)(o) 


(15) 


We can compute contributions to n*^^^(0) from Feynman 
diagrams with topologies shown on Fig. 2 in 

M: 


^(0) “t“ YLxssS ^Zssss 

“h ^Ussss ^^sssss ^Wsff 

d" ^Vf FFFS (16) 

where Hwssss > n^sss v ■ • correspond to topologies (VF, S,...) 
respectively. Structure of these functions is displayed in 
the Appendix. We have regulated IR divergencies of the 
components of (IT^ by adding a mass term to all 
massless sfermionsHI The final result is finite in the limit 
WiR ^ 0- 

We computed all the ingredients in dT51) and expressed 
them via 1-loop integraltH 


A(x) = x(lnx—1), B(x,y) 


1 

J dtln[tx+{l—t)y], (17) 
0 


the 2-loop integrals 


I{x,y,z)=C^ J dfik J d^q 


[P -I- x][g2 -I- y][[k + qY 


-zV 

(18) 


and derivatives of I{x,y,z). C = 167r^p^'^/(27r)‘^, y is the 
regularization scale in DR scheme and e = 4 — 2e. Finally, 
Inx = ln(x/(5^)) where is the renormaliza¬ 

tion scale. Due to supersymmetry, all 1/e^ and 1/e terms 
cancel from 2 so do all diagrams involving counter¬ 
terms. 

Using explicit expressions for /(x, y, z) we find that the 
2-loop contribution to soft mass of (j)i generated by at 
the renormalization scale Q = M is of the following form: 


To hnd 1-loop pole mass of t/i at x = 1 one could solve 
the full propagator equation ml — n(^)(m^) = 0 numeri¬ 
cally. 

The discussed behaviour of m and (HU is not a big 
problem because phenomenology restricts x to the regime 
([R|) which is safe from divergence. 

4. 2-loop self-energy 

To compute 2-loop contributions to soft mass es g ener- 
ated by 0, we apply method developed in [^, 1261 For 
(j)i, we calculate 1-loop self energy n(^)(s) at s = and 
2-loop self energy n(^^(0) at s = 0 and use them to find 


ml.2{M) = KMxffh,2[x), (19) 

where fh, 2 {x) is the complete analytical 2-loop function 
defined in the Appendix. The expansion of fh, 2 {x) for 
X <C 1 is: 

A.2(x) = 1 + F(^2 ^ l)x2 + y^(30^' - 43)x4 + ... (20) 


®The tree-level mass of (f>i is zero. 

^Higgs bosons are assumed to be massless for the simplicity. An 
example of messenger-matter coupling for which Higgs mass term 
is not relevant is Yjj^QD or where ^ are Higgs-type 

messengers. 

^B{x,y) is the standard B{x,y) 1-loop function calculated at s = 

0 . 
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One can check that: (i) fh. 2 {x) is positive for 0 < a: < 1, 
(ii) fh, 2 {x) is divergent at a; = 1, analogously to (nni, and 
(hi) the numerical factor of 3 in (fTOl) agrees with the cor¬ 
responding numerical factor obtained using wave-function 
renormalization method The plot of the function 

fh, 2 {x) is shown on the Fig. [T1 where fh, 2 {x) is compared 
with its approximation (EHl). For X < 0.5 the difference 
between (EHl) and the full result is smaller than 1%. 



Figure 1: The complete 2-loop function fh, 2 i^) defined in (Hsj and 
(I46II (solid line) and its approximation, ll20t . valid for a: <C 1 (dashed 
line). Note that 2 {^) diverges for x 1. 


5. Analytical results: Consistency checks and phe¬ 
nomenological consequences 

To verify our results we have made two basic consis¬ 
tency tests. Namely, we have checked that: (i) in the limit 
X —>• 0 (SUSY limit) the obtained 2-loop contribution to 
soft-masses dni) is zero, and (ii) the sum of 1- and 2-loop 
pole masses m\{Q) = m\ + 2 (Q) does not depend 

on the renormalization scale Q, i.e. 

QdQml{Q) = 0 ( 21 ) 

up to terms of order 0{h^). To this end we have used 
the following 1-loop (3 functions [2^ for the running pa¬ 
rameters: the messenger mass M{Q), F'-term F{Q) and 
superpotential coupling h{Q): 


= hHl, 

(22) 

= -h^F, 

(23) 

= 3h^. 

(24) 

Let us remark that the term (0)n(^)'(0) in 

cial for the cancellation of QdQm\ 

(ITHll is cru- 


In addition, analytical formulas for integrals in 
have been cross-checked using numerical libraries in a wide 


range of (inass)^ parameters: from 10“^^ GeV^ to 10^° GeV^. 
There are a few numerical libraries which can deal with 2- 
loop self energies: BQKASUN 
FIESTA 


SecDec 


CSectors |33l. l34l|. 
and TSIL 371. For our pur¬ 
pose the most appropriate turned out to be the SecDec- 
3.0.6 with Suave routine from the CUBA library (^ . 

We have checked that due to interplay between 1- and 
2-loop contributions to soft masses the apparent domi¬ 
nance of 1-loop contributions is not true for x ~ 1 and it 
turns out that for x ^ 0.9 the total contribution is bounded 
from below, see Fig. [21 instead of having a running direc¬ 
tion caused by the 1-loop contribution. 



X 

CO 

O 



0.80 0.85 0.90 0.95 1.00 
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Figure 2: Comparison of 2-loop (blue, dot-dashed) vs. 1-loop (red, 
dashed) contribution to soft mass for fc ~ 1. The superpotential 
coupling h is fixed to 0.8. The middle curve (solid) shows the sum 
of 1- and 2-loop contributions, ^ 2 ’ g^i^^rated by the 

messenger-matter coupling SJi. 


One can check that the position of the minimum of 
ml 2 on Fig. 12] only mildly depends on the value of h, e.g. 
increasing h from 0.8 to about 1.2 changes its position from 
0.975 to 0.96. 

To show the role which 1-loop function ca and com¬ 
plete 2-loop function (H51) play in the phenomenology of 
the EGMSB models we consider a specific example of o 
in which down-Higgs-type messenger superfield cou¬ 

ples to lepton superfields L, and F: 

W = (25) 

We assume that hij has similar hierarchy to the matrix 
of leptonic Yukawa couplings {ye)ij- One can check by 
direct computation that such superpotential gives extra 
contributions to Ar and to the soft terms at 1- and 2-loop 
level only for sleptons: 

At = -j^3Mxfh,A(^), (26) 
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m 


h,h 


m 


h,TR 


(t^) 4(Mx)V/..2(a^), 


^ (^-ij 2iMx)AAhAx) 

+ (J^^ AMxffhAA- 


(27) 


\ IGTT^y 


(28) 


These soft terms have been implemented iii SARAH- 
4.5.8 [liliE 1^,1^ . which we used to generate SPheno- 


3.3.7 code 5'44l|. The obtained data were further pro¬ 
cessed with the help of SSP-1.2.1 [i^. For simplicity, 
in our analysis we fixed tan/? =10 and chose Mx = 5 x 
lO^GeV such that the lightest Higgs mass m^jO is 123 GeV, 
while masses of gluino g and 1st and 2nd generation squarks 
qi ,2 are about 4 — 5TeV. For h = 0.1, 0.9 and h ~ 2.2 we 
scanned mass spectrum over M in the range 5 x 10® GeV < 
M < 5 X 10® GeV what is equivalent to varying x in the 
range 10“^ < a; < 1. As expected from (l?71) and (1^ 
only 3rd generation slepton masses are sensitive to the h 
coupling. Masses of other particles are similar to those in 
pure GMSB scenario: 


~ — 1.8 TeV, 

4.1 TeV, 


mjro ft! 0.8-1.7TeV, to~± ft 1.5 - 1.7TeV, 


‘‘eh^P-L 


4.4- 5.1 TeV, 
1.7 TeV, 
0.7 TeV. 


s 4.8- 5.1 TeV, 
1.7 TeV, 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 


Dependence of the signed mass of the right-handed stau 
tr i.e. sign(TO|^) on the messenger mass M for 

h = 0.1 and 0.9 is shown in Figs. |3] and |4] respectively. 
When one increases h up to about 2.2 then the depen¬ 
dence on M changes such that is positive for any 
value of M - see Fig. [5] As one can see the interplay be¬ 
tween 1- and 2-loop contributions to soft masses discussed 
in Sec. IV leads to large negative values of when M 
is close to Mq = 5 X 10® GeV. But when M — Mq is of the 
order of 10® — 10®GeV then 2-loop contributions start to 
dominate and is again positive. Such a behaviour is 
apparent especially for h ^ 2.2 (see Fig. [5]). One can see 
that in this case is quite sensitive to the precise value 
of h. The basic difference between discussed scenario and 
the standard GMSB mass pattern is that 3rd generation 
slepton masses can be as well of the order of 500 GeV (as 
in the pure GMSB scheme) as much larger i.e. of the order 
10 — 20 TeV. In addition, contrary to the pure GMSB case, 


^Related SARAH file called EGMSB.m can be downloaded from 
WWW. tjel .us. edu.pl/tools .html. To generate SPheno code please 
use: Start["MSSM"] and then MakeSPheno[lnputFile—>-"EGMSB.ni"] 



Figure 3: Signed mass sign(m^)-y^|m^| of right-handed stau tr as 
a function of the messenger scale M for h = 0.1 and fixed Mx = 
5 X 10® GeV. The subplot displays the behaviour of signed mass for 
M close to Mq = 5 X 10® GeV. The dot-dashed line represents mass 
of tr in pure GMSB scenario (i.e. h = 0). 



Figure 4: Signed mass sign(m^)-\/|m^| of right-handed stau tr as 
a function of the messenger scale M for h = 0.9 and fixed Mx = 
5 X 10® GeV. The subplot displays the behaviour of signed mass for 
M close to Mq = 5 X 10® GeV. The dot-dashed line represents mass 
of Tr in pure GMSB scenario (i.e. h = 0). 


for the fixed value oi Mx slepton masses heavily depend 
on the the value of M. Here they can even be tachyonic 
what could lead to instability of the scalar potential. Fi¬ 
nally, let us note that for M > 10® GeV, what is the most 
common choice of messenger scale, (l?7l) and (1^51) dominate 
over GMSB contributions making masses of 3rd generation 
of sfermions of the order of 10 TeV. Such a hierarchy in 
the slepton spectrum is a simple consequence of choosing 
((Ml) similar to the MSSM leptonic Yukawa couplings. As 
one can see complete 2-loop contributions to soft masses 
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Figure 5: Signed mass sign(m^)-^/|m2| of right-handed stau as a 
function of the messenger scale M for fixed Mx = 5 X 10® GeV and 
h = 2.208, 2.209, ..., 2.219, 2.220 respectively, where the top-most 
line corresponds to h = 2.220. The dot-dashed line represents mass 
of fji in pure GMSB scenario (i.e. h = 0). 


considered here have a large impact on the 3rd generation 
of sleptons. 

6. Summary and Outlook 


suppression by occurs only at 1-loop and not at 2-loop 
level, one might conjecture that na = 0. We also expect 
that the function / 3 (x) behaves similarly to f 2 (x) i.e. it 
is of the order of 1 — 10 for x < 0.9, what in turn could 
suggests that the ratio of 3- to 2-loop contribution is of the 
order of {k^ jk 2 )h? j <0.1. Certainly, from the theo¬ 
retical point of view it would be interesting to obtain even 
approximate 3-loop result and verify those assumptions. 
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Appendix 

Here we collect contributions to self-energy m from 
topologies (W, S',...) defined on Fig. 2 in [2g. We use 
the following notation: z± = M^(l ± x), while c = 
regulates IR divergencies in 2-loop integrals at s = 0. 

The sum of (I55)) - (H5|) is finite in the limit c —>■ 0. The 
components of (HU) are expressed in terms of functions 
{Wssss,^SS,Yssss, ■ ■ ■) defined in the Sec. IV.A and 
IV.B of 0. 


In this letter, we have presented full 2-loop calcula¬ 
tion in the EGMSB model with one coupling between 
messenger and matter. We have shown that full 2-loop 
analytical contributions to soft masses are important in 
the regime x = F/M"^ ~ 1, which is realized e.g. in the 
Low-Scale Gauge Mediation models with messenger mass 
M ~ 10® — 10®GeV. It turned out that for x ~ 1 the 
model exhibits novel behaviour with respect to the stan¬ 
dard GMSB case. It has been shown on a specific example 
that considered here complete 2-loop effects can signifi¬ 
cantly affect slepton spectrum, what opens up new sce¬ 
narios beyond the standard GMSB scheme. 

Let us shortly comment on the remaining contributions 
from higher loops. Taking into account suppression of 1- 
loop contribution with respect to 2-loop contribution, the 
soft mass can be written in the following form: 


m 


2 

h 


(Mx)2 


+ (y^) ksx'^^ f3{x) + ..., (35) 


where ki = —1/6, k 2 = 3, k^ is unknown numerical factor 
of 3-loop contribution, while fi{x) are functions normal¬ 
ized such that /i(0) = 1. Finally na > 0 is unknown power 
of X which appears in the 3-loop contribution. We are not 
aware of any strong argument that na = 0 and /a (a;) is of 
the order 1 for x < 0.9, but from the fact that in the 


I^^^^Wssss 2IT5555(2:-|-, Z—, C, c) 

+2IEssss(c, c, c,z-) + Wssss{z-,z-,c, c) 
+2ITssss(c, c, c, z+) + H/ssss(z+, z+, c, c), (36) 


+Xsss(,z+, Z-,c) + Xsssic, c, Z-) 

+Xsss(,z-,z-,c) + Xsss{c, c, z+) 

+-^sss(^-i-, z+,c) + 2Xsssic, c, c), (37) 


^4^^Yssss Z— , Z-^-, c) 

-2Xssss(c, z+, z-,c)+ Yssss{z-,c, c, z+) 
+Yssss{z+,c,c, Z-) 

+2Yssss{c, Z-, z-, c) + 2Yssss{z-,c, c, c) 
+Yssss{z-,c, c, z_) + 2Yssss{c, z+,z+,c) 
+2issss(2:+, c, c, c) + Yssss{z+,c, c, 2 :+), (38) 
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^^^^Zssss ^ssss{^: Z—^ 

-Zssss{c, Z-, z+,c) - Zssss{c, z+,c, Z-) 

-Zssss{c, Z+,Z-,c) - Zssss{z-,C, C, Z+) 
-Zssssiz-,C, Z+, c) - Zssss{z+, C, C, Z-) 
-Zssssiz+,C, Z-,c) + Zssss{c, Z-,C, Z_) 

+ Zssss{c, Z-,Z-,c) + Zssss{z-,C, C, Z-) 
+Zssss{z-,c, z-,c) + Zssss(c, z+, C, Z+) 

+ Zssss(c, z+, Z+, c) + Zssss(z+,C, c, z+) 
+Zssss(z+,c,z+,c), (39) 

4 

J^^Ssss ‘^SSs(c, Z—, Z-f-) 

+Ssss(c, Z+, Z-) + Ssss(c, Z-,Z-) 

+Ssss{c, z+, z+) + 4Ssss(c, c, c), (40) 

4 

^4^nf/gggg UssSsi,Z—: Zj^') 

-Ussss{z-,C, Z+,c) - Ussss{z+,C,C, Z-) 
-Ussss{z+,c, z-,c) + 4Ussss{c, z-,c, c) 
+Ussss{z-,c, c, Z-) + Ussss{z-,c, z_,c) 
+4Ussss(c, z+,c, c) + Ussssiz+,c, c, z+) 
+Ussss{z+,c,z+,c), (41) 

4 

^4^2 ^SSSSsi^^ Z—^ Z-\-^ C, c) 

+ Vsssss{c, Z+,Z_,C, c) + Vsssss{z-,C, c, c, z+) 

+ Vsssss(^+,C, C, C,Z-) + Vsssss{c, Z-,Z-,C,c) 

+ Vsssss(z-,C, C,C, Z-) + Vsssss{c, z+,z+,c, c) 
+Vsssss(^+,c, c, c, 2 :+), (42) 

4 

■j^4Jwsff = 4Wsff{c, c, 0, z) 

+Wsff{z-, z+, 0,0) + Wsff{z+, z-,0, 0) 
+Wsff{z-, Z-,0, 0) + Wsff{z+, 2+, 0,0), (43) 

4 

'J^4~^^VsssFF s SF f{,Z— C, C, 0, z) 

+ 2VssSFf{z+, C, C, 0, z) - VssSFFic, Z-,Z+, 0,0) 
-Vsssff{c, z+,2:_,0, 0) + Vsssff{c, z-,z-,0,0) 
+Vsssf_f(c, 2:+,2:+,0,0), (44) 

2 

■j^41vffffs = 4Vffffs{0,z,z,0,c) 

+2Vffffs{z, 0,0, z, c) + Vffffs{z, 0,0,0, z_) 
+Vffffs( 2, 0,0,0, z+). (45) 


The complete analytical form of the function fh, 2 {x) 
which enters the 2-loop contribution to soft mass (HU) is: 
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fh,2{x) = 

1 


12{x + 1)^ (a;^ — 2x -I- 2) a:^Li 2 (l — x) 


72x2(1 _ a;2)2 _ 

-12(x - 1)^ (x^ -h 2x -h 2) x^Lia 

—24 (x^ -I- l) x^Li 2 ^ in(2x)artanh(x) 

— 12x^ inxin(x + 1) — 24(x^ — l)(x — l)x^C(2) 

-|-48x^ in(2)artanii(x) -I- 24(x^ — l)x^ 

—6(x^ — 1) (l9x^ — 20) in (l — x^) 

— 12x(x^ — 1) (7x^ — lO) artanh(x) 

-1-12 (x^ -I- 2) x^ inxin (l — x^) 
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6x2(1 — x2) 
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tionai materiai in www.tjel .us . edu.pl/tools .html. 
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